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We study numerically the dynamo transition of an incompressible electrically conducting fluid 
filling the gap between two concentric spheres. In a first series of simulations, the fluid is driven 
by the rotation of a smooth inner sphere through no-slip boundary conditions, whereas the outer 
sphere is stationary. In a second series a volume force intended to simulate a rough surface drives 
the fluid next to the inner sphere within a layer of thickness one tenth of the gap width. We 
investigate the effect of the boundary layer thickness on the dynamo threshold in the turbulent 
regime. The simulations show that the boundary forcing simulating the rough surface lowers the 
necessary rotation rate, which may help to improve spherical dynamo experiments. 
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I. INTRODUCTION 

Dynamo theory describes the generation of magnetic 
fields in flows of conducting fluids, for example in stars or 
planetary interiors. Several experiments have been built 
in order to reproduce a dynamo in the laboratory. Af- 
ter the first experiments in Riga [l| and Karlsruhe Q, 
more recent experiments 0, 13 implement flows less con- 
fined to specific shapes than the first experiments, so that 
the effect of turbulence on dynamo action is more rele- 
vant. Nevertheless, the main contribution to magnetic 
field generation in these experiments is still assumed to 
come from the time averaged part of the flows. The nu- 
merical investigations presented in this paper are directly 
applicable to spherical experiments like those in Mary- 
land ||. 

The spherical Couette flow with a stationary outer 
sphere is mainly a differential rotation driven by the in- 
ner sphere, together with a jet in the equatorial plane 
in which fluid is centrifuged radially outward. A flow 
along the rotation axis connecting the poles of the two 
spheres brings the fluid back to the inner sphere. The 
flow thus consists of differential rotation and two loops in 
the meridional circulation. It is helical with opposite he- 
licities in the two hemispheres and is topologically equiv- 
alent to the s2tl flows studied by Dudley and James [5| 
which are known to be dynamos. The spherical Couette 
flow on the contrary does not generate magnetic fields 
for Reynolds numbers Re up to a critical value where 
the flow becomes unstable [6]. It therefore seems possi- 
ble that turbulent spherical Couette flow is a small scale 
dynamo. In this type of dynamos, at magnetic Prandtl 
numbers Pm with Pm ^> 1, the largest production of the 
magnetic field occurs at the resistive scale and the critical 
magnetic Reynolds number Rm increases with increas- 
ing hydrodynamic Reynolds number. For Pm <C 1, the 
largest production of magnetic field occurs at some scale 
larger than the viscous scale @, 0, so that the critical 
Rm is independent of Re for Re large enough, because 
increasing Re then adds small vortices with magnetic 



Reynolds numbers less than 1. The magnetic field shows 
nothing but diffusive dynamics on these length scales so 
that they do not contribute to either creation or mixing 
of magnetic field. Ref. 6 found Rm c to increase with Re 
for the entire investigated parameter range, which further 
adds to the suspicion that spherical Couette flows could 
be small scale dynamos. Since the Reynolds number Re 
in numerical simulations is limited, the magnetic Prandtl 
number Pm has to be adjusted to values far above the 
Pm of liquid sodium (which is the liquid commonly used 
in experiments) in order to achieve high enough magnetic 
Reynolds numbers. There is always a need to extrapolate 
from numerical results to larger Re, so that it is always 
important to understand the dependence of the dynamo 
threshold on Re and Pm [7H9j. 

An experiment is under construction which realizes 
spherical Couette flow in liquid metals [lOj. One mo- 
tivation for the present study is to predict whether this 
experiment will be able to sustain a self-generated mag- 
netic field. This question will be investigated in section 
IIII1 Because of the pessimistic answer obtained in this 
section, section [TV] considers the possible improvement 
obtained by welding blades on the inner sphere in or- 
der to strengthen the coupling between the fluid and the 
rotation of the inner sphere. 



II. THE MATHEMATICAL MODEL 

The system under investigation consists of two concen- 
tric speres with radii Ri and R Q forming a gap of width 
Ri — R = d and aspect ratio Ri/R = 1/3. The fluid in 
the gap is driven by the rotation of the inner core, which 
rotates at angular frequency Qi about the z— axis, while 
the outer boundary is at rest. The fluid is characterized 
by its kinematic viscosity v and its magnetic diffusivity A. 
The evolution of the magnetic field in an incompressible 
electrically conducting fluid is described by the induction 
equation. The system is governed by two numbers, the 
Reynolds number Re, and the magnetic Prandtl number 
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Pm or alternatively the magnetic Reynolds number Rm 
defined by 



Re 



Pm 



v 
A' 



Rm = Re Pm. 



(i) 



From now on, only adimensional variables will be used 
in this paper. Choosing as units of time and length the 
reciprocal of the inner core's rotation rate and the gap 
width, the induction equation reads for non-dimensional 
magnetic and velocity fields B and v: 



d t B + Vx(Sxv) 



1 



Re Pm 



V 2 B , V-B = 0. (2) 



The velocity field itself is a solenoidal vector field deter- 
mined by the Navier- Stokes equation and the continuity 
equation: 



d t v + (v • V)v 



-V$ + — V 2 v 
Re 



V-u = 0. (3) 



Eq. (|3|) contains a volume force F and $ stands for 
the pressure variable. Eqs. (j2j) and (j3]) describe the kine- 
matic dynamo problem, in which one assumes the system 
to be near the onset of magnetic field generation, so that 
the magnetic field strength is small and the Lorentz force 
is negligible in equation (j3j). 

The parameters Re and Rm depend on the inner 
core's rotation rate. More revealing parameters are the 
Reynolds numbers Re and Rm based on the rms velocity 
defined as 



= ( / \v 2 dV) 



and 



Re = Re v r 



Rm = Re Pm 



(4) 



(5) 



where V is the volume of the shell, brackets denote time 
average and E^ n the kinetic energy. 

In the following the effect of different surface proper- 
ties on the dynamo threshold are compared. The inner 
and outer spheres have radii vi = 1/2 and r Q = 3/2, 
respectively. The inner sphere and the space surround- 
ing the outer sphere are assumed to be insulating. Two 
different types of boundary forcing are implemented in 
order to simulate both smooth and rough surfaces. For 
the smooth surface no-slip boundary conditions 
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are chosen together with F — in eq. (j3]). In the second 
case, the inner boundary is assumed free slip and the 
fluid is driven in one tenth of the gap width by a toroidal 
force field F of spherical harmonic degree I = 1 and order 
m = 0, given in spherical polar coordinates (r, 0, cp) by: 



F= 1 - 
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(7) 



This forcing qualitatively reproduces the flow driven 
by an inner sphere with blades of height one tenth the 
gap size mounted along meridians. The non-dimensional 
rotation rate of the inner sphere, has to be found a 
posteriori in these simulations by computing 
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dS sin 2 6(v lp (r = r i ,6,<p,t)) (8) 



with which a Reynolds number Re', analogous to Re for 
the no slip boundaries, is defined as 



Re' = Re Q' 



(9) 



In the following Re is going to be increased up to 
1.67 • 10 4 for the no-slip boundary conditions and up to 
2.5 • 10 3 in the simulations utilizing the volume force. 
The magnetic Prandtl number is of order unity for the 
simulations just at the onset of magnetic amplification. 
The numerical method is a spectral method, in which the 
fields are expanded in spherical harmonics and Cheby- 
chev polynomials [11] . 



III. NO-SLIP 

A. Hydrodynamic Characteristics 

The hydrodynamic properties of the spherical Cou- 
ette flow with stationary outer sphere have already been 
investigated in detail in [f| [l2|, [l3|. At low Reynolds 
numbers the basic spherical Couette flow is axisymmet- 
ric and we find a critical Reynolds number of Reh = 1500 
(Reh = 105) beyond which small non-axisymmetric per- 
turbations increase and an instability develops as a prop- 
agating wave on the equatorial jet with a dominant az- 
imuthal wavenumber m = 2, which is in agreement with 
[|, [l3|. At Re s = 2800 (Re s = 178) a second transi- 
tion occurs. Beyond this value, the amplitudes of odd 
wavenumbers m of the kinetic energy develop as well 
and the power spectrum begins to flatten and approaches 
power laws in m _5//3 and / -5 / 3 , indicative of Kolmogorov 
turbulence (see Figure [1]). Large scales of the velocity 
field, however, have comparable spectra for all Re. The 
bottom panel shows the power spectra of an estimate of 
the turbulent rate of strain of spherical harmonic order 
1, which is simply the kinetic energy multiplied by l 2 and 
has its maximum at the viscous scale. At Re = 970 the 
inertial range reaches up to i« 40. 

The wavenumber of the most unstable mode depends 
on the aspect ratio of the shell and switches from m — 2 
to m = 3 at an aspect ratio close to the one chosen here 
pjl . It is possible to generate flows with m = 3 as domi- 
nating mode by starting from an equilibrated solution at 
Re > Re s and lowering Re to some value in between Reh 
and Re s . 

Figure [2] (top) shows the thickness of the boundary 
layer of near the inner core. In order to produce this 
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FIG. 1: (Color online) Power spectrum of kinetic energy plot- 
ted against spherical harmonic degree m (top) and the turbu- 
lent rate of strain versus spherical harmonic order I (bottom) 
for Re = 300 (black continuous line) and 970 (thick red dashed 
line), together with the Kolmogorov power laws m~ 5 ^ 3 and 
I 1 / 3 (blue dashed line). 



figure, Up has been averaged over spherical surfaces and 
the boundary layer thickness denned as the distance from 
the inner sphere at which this averaged velocity drops to 
the arithmetic mean of its value at the boundary and its 
value averaged over radius. At high Re the evolution of 

i/2 

the boundary layer thickness is approaching Re 1 , as 
one would expect from theory. 

Figure [2] (middle) shows the dependence of the dimen- 
sionless kinetic energy on Re, which converges to a con- 
stant in the limit of Re -1 — >• 0. The viscosity becomes 
irrelevant for high Re and the only remaining control pa- 
rameters entering the dimensional kinetic energy are the 
density and so that the kinetic energy scales with fi? 
and the dimensionless kinetic energy becomes constant. 

A quantity of direct relevance to experiments is the 
energy dissipated in the flow. The energy budget, ob- 
tained by taking the scalar product of eq. (|3j) with v and 
integrating over space, reads 



d t J ^v 2 dV = ^r-^ J (d iVj ) 2 dV + J F-vdV (10) 
with 

T = -j o d0rS Sin2 (?rV<p ~ y) , (11) 

evaluated on the inner boundary r = T{ , being the torque 
on the inner boundary. F = for the simulations in this 
section so that the time averaged torque (r) is directly 
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FIG. 2: (Color online) Boundary layer thickness ds plotted 
against Re -1 / 2 (top), kinetic energy plotted against Re -1 
(middle), and energy dissipation rate plotted against Re (bot- 
tom). The vertical dashed lines mark Reh and Re s . 



related to the energy dissipation e by e = (r)/Re. As 
shown in fig. [2j the dissipation scales as Re -0,62 . The 
bracket in the integrand in eq. (fTTj) is dominated by the 
derivative d r v^ which can be estimated as v^qi /(Ib, where 
ds is the boundary layer thickness and v re \ the velocity of 
the inner boundary relative to the bulk fluid. Since ds oc 
Re -1 / 2 , we have e ex Re~ 1 ^ 2 v re \. If the motion of the 
inner core was completely decoupled from the rotation 
of the inner boundary, v re \ would be independent of Re 
and e ex Re 1 . There is of course some entrainment of 
the fluid by the rotation rate of the inner sphere and one 
finds an exponent of -0.62 instead of -0.5. 



B. Dynamo transition 

The dynamo transition for the spherical Couette flow 
with stationary outer sphere has already been computed 
in [6[ with different magnetic boundary conditions and 
an aspect ratio of r^/r = 0.35. These results will be 
compared with ours in the last section. Figure [3] shows 
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the dynamo simulations in the (Re, Rm)— plane. Aster- 
isks mark working dynamos and dots are failed dynamos. 
Linear interpolation between the growth rates computed 
at those points allows us to find the locus of zero growth 
rate. The onset of magnetic field amplification obtained 
in this way is indicated by the dashed line. The mag- 
netic fields are dominated by modes with an azimuthal 
wavenumber of m = 2. As mentioned above, suitable 
initial conditions lead to magnetic fields dominated by 
modes with m = 3. The onset for these dynamos is in- 
dicated by the thick solid line. The two lines coincide 
for Re > Re s . The straight line indicates Pm = 1. Ex- 
cept for the highest Re the magnetic Prandtl number 
is always larger than one. In agreement with [(| there 
are no working dynamos for axisymmetric flows. For 
Reh < Re < Re s the critical magnetic Reynolds num- 
ber Rm c for m = 2 first increases, reaches a maximum 
and finally decreases, until at the second transition odd 
modes of the kinetic energy become unstable and the dy- 
namo threshold again increases. 
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FIG. 3: (Color online) Dynamo transition for magnetic fields 
dominated by modes with m — 2 (black dashed line) and 
m = 3 (blue bold line). Failed dynamos are indicated by red 
dots, working dynamos by blue asterisks. The vertical dashed 
lines seperate the different hydrodynamic regimes: Axissym- 
metric flow (left), first nonaxisymmetric instability (middle) 
and turbulent regime (right). 



of Re, but neither the frequency of the oscillation of the 
magnetic energy nor the growth rate agree exactly. How- 
ever, the amplitude of the m = 4 contributions to the 
kinetic energy are only one order of magnitude smaller 
than the dominant mode and they introduce a slight time 
dependence in the codrifting frame of reference, so that 
the simulations of the drifting snapshot can only catch 
the qualitative features of the full simulation. 
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FIG. 4: (Color online) Power spectrum of the kinetic energy 
(black continuous line) and the magnetic energy (red dashed 
line) for Re = 130 (top) and Re = 160 (bottom) 



The non-monotonous dependence of Rm on Re in the 
interval Reh < Re < Re s appears to be related to two 
other features of the dynamo: The drift frequency of the 
unstable mode and the dominating wavenumber in the 
magnetic field. 

The velocity field for Reh < Re < Re s is dominated by 
one azimuthal wavenumber and its harmonics. Time se- 
ries of the radial velocity field at a fixed point show that 
the phase velocity of the propagting wave diminishes for 
increasing Re. Similarly to previous studies of kinematic 
dynamos [3 EH , it seemed useful to run simulations in 
which the velocity field is a snapshot taken from the full 
simulation, and which is set to drift at arbitrary phase 
velocities in order to investigate the dependence of the 
growth rate on the phase velocity of the instability. It 
turned out that the critical magnetic Reynolds number 
increases with the phase velocity, as does Rm c in Figure 
[U The time evolution of the magnetic energy using the 
artificially drifting velocity field shows the same charac- 
teristic superposition of an exponential growth and an 
oscillation as in the dynamic simulations in that region 
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FIG. 5: (Color online) Power spetrum of the kinetic energy 
(black solid line) and the magnetic energy (red dashed line) 
for Re = 130 in a simulation in which modes with m — 3 
dominate. 

In Figure [H power spectra of the kinetic and magnetic 
energy of the simulations at Re = 130 and Re = 160 are 
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shown. In the first case only even wavenumbers of the 
magnetic energy are amplified, whereas odd wavenum- 
bers decreases in time. The opposite is seen in the sec- 
ond case, where odd wavenumbers of the magnetic en- 
ergy are amplified and even wavenumbers decrease. The 
different symmetry of the magnetic field seems to be re- 
sponsible for the decrease of the dynamo threshold, since 
this decrease correlates with the appearence of the odd 
wavenumbers in the magnetic field. In addition, Figure 
[5] shows the spectra of the simulations with m — 3 as the 
dominant wave number, which have for Reh < Re < Re s 
the same symmetry and accordingly the dynamo thresh- 
old is monotonously increasing. 

For Re > Re s , the critical Rm first increases but then 
reaches a plateau at Rm = 600. Pm is less than 3 
throughout the plateau region. It will now be argued 
that even more turbulent flows will not alter the critical 
Rm significantly. 
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FIG. 6: (Color online) Power spectrum of magnetic energy of 
spherical harmonic degree m (top) and order / (bottom) for 
Re = 300 (black continuous line) and 970 (red dashed line). 

In the plateau region, the Reynolds number is high 
enough for the kinetic energy spectrum to decay as a 
power law as a function of both m and I with an exponent 
of -5/3. The phenomenolgy developed for small scale dy- 
namos in Kolmogorov turbulence [7|, [8[ therefore applies. 
According to this phenomenology, Rm can increase as a 
function of Re if Pm is larger than 1 and the eddies near 
the viscous length scale generate the magnetic field. If 
either Pm is small or the field generation occurs whithin 
the inert ial range or at the integral scale, Rm is indepen- 
dent of Re. In spherical Couette flow, the large scales are 
already capable of dynamo action. Not surprisingly, the 
spectra of magnetic energy peak at small m and I (see 
Fig. [6j). Only the computation at the largest Re raises 



doubts because the spectrum of magnetic energy reaches 
a second local maximum at I = 11 in this case, which 
nevertheless is still within the inertial range (see Fig. [T]) 
and at any rate below the global maximum at I = 1. Pm 
is already slightly below 1 in this computation so that a 
further increase in Re should only add a high wavenum- 
ber tail to the kinetic energy spectrum corresponding to 
eddies with too small a magnetic Reynolds number to 
contribute to the dynamo effect. In order to further test 
whether the small scales are responsible for magnetic field 
generation, eqs. (03]) were solved simultaneously, but the 
axisymmetric components of the velocity were removed 
before the induction term in eq. (j2j) was computed. The 
magnetic field decayed in this simulation, which shows 
that the turbulent eddies in this velocity field are unable 
to support the magnetic field by themselves. 

Since the magnetic field is generated at large scales 
even for Pm « 1, an increase in Re and concomittant 
decrease in Pm adds scales which contribute neither to 
the creation nor the destruction of the magnetic field. 
It is concluded that Rm = 600 is the critical value of 
Rm in more turbulent flows as well. In an experiment 
with Pm = 10 -5 this corresponds to Re = 6 x 10 7 . 
For Reynolds numbers this large, i^kin is nearly inde- 
pendent of the Reynolds number according to fig. [2] and 
equal to 0.023, from which one deduces v rms = 0.058 and 
Re ~ 10 9 . For the sodium experiment in Maryland with 
a gapwidth of d — lm and a viscosity v w 10 -6 m 2 /s, 
this corresponds to a rotation period of about only 6ms. 



To conclude this section, we compare the stability di- 
agram in fig. [3] with the results of similar studies. All 
numerical simulations report a general increase of the 
critical Rm with the Reynolds number for Pm larger than 
1, an effect which is also found in the analytical model 
of ref. [16[ . The dependence of the dynamo threshold on 
the Reynolds number is not monotonous, however, and 
other published results have in common with fig. [3] a 
maximum in the critical Rm for Pm around 1. This fea- 
ture is seen in refs. [8, 9]. Explanations for this effect 
have been attempted in refs. [l?], G3 who attributed it 
to either the bottleneck effect in turbulent spectra or the 
helicity at the viscous scale in the velocity field. 

Despite the similarities, it is doubtful whether there is 
a universal mechanism for the dependence of the dynamo 
threshold. For example, the plateau in Rm is reached for 
values of Pm ranging from 0.1 to 1 across the different 
studies, so that there are variations by an order of mag- 
nitude in the Pm at which the critical Rm is maximum. 
There is also no indication of a bottleneck in the spec- 
tra presented here which shows that the bottleneck is 
not necessary for the observed variation of the dynamo 
threshold. It remains an open question why there is a 
maximum in the critical magnetic Reynolds number as 
shown in fig. [3] in the spherical Couette system. 
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IV. ROUGH SURFACE 
A. Hydrodynamic characteristics 

The critical Re found in the previous section is chal- 
lenging to realize in experiments. The problem arises 
from the fact that the inner sphere needs to rotate faster 
in order to increase Re. But a faster rotation rate also 
thins the boundary layer at the inner sphere and reduces 
the volume of fluid directly coupled to the boundary mo- 
tion. An obvious remedy is to mount blades on the inner 
surface. The boundary layer cannot be smaller than the 
surface roughness which ensures a better coupling of the 
fluid to the inner sphere at high rotation rates. In the 
following, the blades are modelled by a volume force as 
described in section ITT1 

Due to the increased momentum transport, turbulence 
develops at lower rotation rates. The first nonaxisym- 
metric instability arises at Reh = 425 (Reh = 95) with a 
dominant wavenumber of m = 2, and beyond Re s = 465 
(Re s = 108) odd wavenumbers increase, too. The evolu- 
tion of the boundary layer thickness is shown in Figure [71 
(top). As expected, the thickness of the boundary layer 
does not drop below one tenth of the gap width and ap- 
proaches the thickness of the forced layer at large Re. 
Since the thickness of the layer in which momentum is 
injected into the fluid remains nearly constant, indepen- 
dend from the rotation rate, the kinetic energy increases 

1/2 

with Re. It does so approximately in Re , which is 
shown in Figure [71 (middle). 

The torque on the inner sphere is zero for free slip 
boundary conditions so that one deduces from the energy 
budget ([10)) that the dissipation rate is given by e = 
J F ■ vdV. As can be seen in the bottom panel of fig. [7J 
e behaves differently depending on whether the boundary 
layer is thicker than the forced layer or not. The interval 
of Reynolds number in each regime is too small to deduce 
a law relating e with Re'. 
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FIG. 7: (Color online) Boundary layer thickness du plotted 
against Re _1//2 (top), kinetic energy plotted against Re 
(middle), and energy dissipation rate plotted against Re' (bot- 
tom). The vertical dashed lines show transitions to different 
regimes. Results obtained for smooth and rough surfaces are 
indicated by the symbols x and +, respectively. 



B. Dynamo transition 



a- = 9-4. 



The diagram of the dynamo threshold plotted in the 
(Re, Rm)— plane in Figure [8] is similar to Figured! There 
is no dynamo action for axisymmetric flows. Similar to 
the simulations with no-slip conditions, even wavenum- 
bers (m = 2 and harmonics) increase in the magnetic 
energy spectra, whereas the total magnetic energy is in- 
creasing exponentially with superposed oscillations. In 
contrast to the simulations of the previous section, it was 
not possible to find a case in which the odd wavenumbers 
dominate the magnetic spectrum. Like in the simulations 
with no-slip conditions the dynamo threshold shows a 
plateau at Rm = 600. In experiments with Pm = 10 -5 
this corresponds to Re = 6 x 10 7 . By extrapolating the 
kinetic energy in Figure [7] (bottom) up to this Re, one 
gets £?kin ~ 300 and « 6.6, or Re' ^ 9 x 10 7 with 




FIG. 8: (Color online) Dynamo transition (black dashed line) 
with the same symbols as in figure [3] 



7 



\ \ 

! ! 

\ i 




y\v^ ****** 

/ 

/ 






— smooth 






— rough 






C.Guervilly 






P. Cardin 



10 3 10 4 

Re 



FIG. 9: (Color online) Dynamo threshold of both surface 
types and results of 



V. CONCLUSION 

Laminar spherical Couette flow with a stationary outer 
sphere is axisymmetric and consists of a differential ro- 
tation superimposed on a meridional circulation, which 
flows from the inner to the outer sphere in the equato- 
rial plane, continues towards the poles of the outer sphere 
along the outer boundary in each hemisphere, and finally 
returns to the inner sphere. Dudley and James [5[ showed 
that flows of this type are capable of dynamo action. At 
the magnetic Reynolds numbers accessible by numerical 
computation, laminar spherical Couette flow nonetheless 
is not a dynamo due to the unfavourable velocity pro- 
file. For instance, the outward flow is concentrated in 
a narrow jet in the equatorial plane. Unstable spherical 
Couette flow on the other hand readily generates mag- 
netic field. However, the field generation process still 
occurs at large scales, the small turbulent scales on their 
own do not support the magnetic field in these dynamos. 

The numerical data allow us to extrapolate to more 
turbulent flows than the simulated flows if the data are 



plotted in terms of Reynolds numbers, magnetic and hy- 
drodynamic, based on the rms velocity rather than the 
boundary velocity. The data are more difficult to inter- 
pret and extrapolate if they are given as functions of the 
Reynolds numbers computed with the inner boundary 
velocity as seen in Figure [9] This plot also includes the 
data from ref. [6[ , which should be compared to our sim- 
ulations with smooth boundaries. The curves are broadly 
similar apart from a shift by roughly a factor 1.5 along 
the Re— axis, which must be attributed to the different 
aspectio ratio and the different magnetic boundary con- 
ditions. 

If the flow is driven by the volume force which simu- 
lates blades of height one tenth of the gap size mounted 
on the inner sphere, the dynamo threshold expressed as 
critical magnetic Reynolds number based on the rms ve- 
locity is reduced by one sixth. Such a change is plau- 
sible because the topology of the flow is the same as 
for smooth boundaries, but the kinetic energy is more 
evenly distributed in space. Both the boundary layer 
and the equatorial jet are thicker. From an experimen- 
tal point of view, the much more important improve- 
ment brought about by the blades is the reduction of the 
rotation rate of the inner sphere at the dynamo onset. 
The critical Reynolds number based on the rotation fre- 
quency of the inner sphere is reduced by a factor 10 by 
the blades. Using the results obtained with the simu- 
lated blades and the numbers for the Maryland experi- 
ment with a gapwidth of one meter and the viscosity of 
liquid sodium v « 10 -6 ra 2 /s, a Reynolds number num- 
ber of Re = 9 x 10 7 is reached for a rotation period of 
the inner sphere of 0.07 s. A spherical Couette dynamo 
is more readily attainable with blades mounted on the 
inner sphere because they improve the coupling between 
the fluid and the motion of the inner sphere and at the 
same time modify little the structure of the flow. 
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